Graphene spintronics: the role of ferromagnetic electrodes 
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We report a first principles study of spin-transport under finite bias through a graphene- 
ferromagnet (FM) interface, where FM=Co(ffi), Ni(lll). The use of Co and Ni electrodes achieves 
spin efficiencies reaching 80% and 60%, respectively. This large spin filtering results from the mate- 
rials specific interaction between graphene and the FM which destroys the linear dispersion relation 
of the graphene bands and leads to an opening of spin-dependent energy gaps of ~ 0.4-0.5 eV at the 
K points. The minority spin band gap resides higher in energy than the majority spin band gap 
located near Ef, a feature that results in large minority spin dominated currents. 



The field of spintronics, or magneto-electronics, utilizes 
the spin degree of freedom of electrons and their inherent 
magnetic moment to influence or control the properties of 
a circuit. Within this field much effort has focused on de- 
veloping interfaces, commonly a non-magnetic metal or 
insulator in contact with a ferromagnet (FM), that ex- 
hibits a large spin-polarized interface resistance £Q. Ide- 
ally, such a spin filter would allow electrons of only a 
single spin component to conduct. 

Graphene, a 2D lattice of C atoms, is a gapless mate- 
rial with linear dispersion electronic bands joining at the 
Fermi level (Ep) in conical (Dirac) points located at the 
K points in the Brillouin zone (BZ) [2J. It has received 
much attention due to it's exceptional properties [3J, in- 
cluding zero effective mass carriers with extremely large 
mobilities, and is poised to play a role in the future of 
nanotechnology. Among other qualities, graphene has 
weak spin-orbit interaction due to the low atomic number 
of C resulting in long spin-flip scattering lengths. Hence, 
graphene is a promising material for applications in spin- 
tronic devices, where one can exploit graphene's unique 
electronic properties within the context of magneto- 
electronics. 

Generating and injecting a spin-polarized current into 
graphene is of vital importance to the development of 
graphene-based spintronics. Graphene nanoribbons, un- 
like pure graphene, are theoretically predicted to possess 
a local magnetic moment at the zigzag edges[3], but a ma- 
jor limitation arises in the difficulty of reliably fabricating 
well defined low-disorder edges. Thus, efficient spin in- 
jection into graphene is required for the realization of a 
prototypical spintronic device. A previous first principles 
study showed extremely large spin filtering efficiencies 
for FM|graphene(Gr)|FM junctions (FM=Co, Ni)|5], in 
which the current flow was oriented perpendicular to the 
graphene. In this way, the spin-polarized current is pri- 
marily dominated by the inter-layer coupling, namely van 
der Waals interactions, between graphene sheets, rather 
than the characteristic graphene states. Moreover, the 
current-in-plane geometry, with transport occurring par- 
allel to the graphene, is the most common experimental 
and theoretical device architecture 6-9 . In such sys- 



tems the source and drain electrodes are comprised of 
FM-covered graphene (because metallic contacts are de- 
posited on top of graphene) which, depending on the 
nature of the chemical bonding, can hybridize and re- 
sult in a complex electronic structure. It was previously 
shown that graphene placed in contact with Co or Ni 
strongly hybridizes leading to significant modifications 
of the graphene bands [5]. Demonstrated by a similar 
systemflQ^, it can be inferred that the interaction be- 
tween graphene and a FM is very sensitive to the par- 
ticular atomic configuration at the interface of the two 
materials. Hence, it is crucial to properly characterize 
the detailed atomic structure, in order to obtain the cor- 
rect electronic states at the interface and accurately ana- 
lyze the spin-polarized transport properties of the device. 
Given the nature of this problem, one must employ atom- 
istic ab initio modeling for an accurate treatment of the 
chemical interaction at the contact. 

A nano-structure was constructed to model the inter- 
face between a source or drain electrode (i.e., the FM- 
covered graphene) and a pure graphene channel, as shown 
in[TJ This interface will be hereafter referred to as Gr|FM. 
This particular choice in system geometry was motivated 
by calculating interface properties that are independent 
of the device length. Thus, whether considering a very 
long (diffusive transport) or very short (ballistic trans- 
port) graphene channel, we expect the results presented 
in this work to remain valid at the interface. We con- 
sider both Co(lll) and Ni(lll) as the FM in contact 
with graphene, forming our electrodes. For graphene- 
based spintronics, Co(lll) and Ni(lll) are excellent can- 
didates for FM contacts since their in-plane lattice con- 
stants nearly perfectly match that of graphene[5], with 
experimental mismatch values of 1.8% (Co) and 1.3% 
(Ni). 

In this Letter, first principles density functional the- 
ory (DFT) total energy and non-equilibrium transport 
calculations were carried out to study the spin-polarized 
electronic structure and spin-dependent transport prop- 
erties of Gr|FM interfaces under finite bias. In particu- 
lar, the atomic structure of the FM-covered graphene was 
fully relaxed by DFT total energy calculations. Given the 
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FIG. 1: (a) Diagram of the Gr|Co(lll) and Gr|Ni(lll) inter- 
face. The dotted lines indicate that the left and right leads, 
extending to ±00, consist of FM-sandwiched graphene (left 
lead) and pure graphene (right lead). The top-site and hollow- 
site C atoms forming the graphene are shown in blue and 
green respectively. Periodic boundary conditions are assumed 
in the plane perpendicular to current, i.e., the system extends 
infinitely along the x- and z-directions. (b) View of the left 
lead in the x-y plane showing the graphene sitting on the 
FM. The top-site C is directly above the FM atom, while the 
hollow-site C is located at the hollow site, (c) View of the left 
lead in the y-z plane indicating the optimized graphene-FM 
distance do- The dotted lines shown in (b) and (c) delimit the 
supercell box used for the electronic structure calculations. 



optimized atomic coordinates, the spin-dependent band 
structure was analyzed. It was found that the strong 
hybridization between graphene and the FM destroys 
the linear dispersion relation of the graphene bands and 
opens spin-dependent band gaps at the K point (similar 
to previous work [S])- This FM- induced band gap open- 
ing of the graphene states results in spin-polarized cur- 
rents that are minority (MIN) spin dominated due to the 
majority (MA J) spin band gap residing near Ep. This 
electronic feature leads to spin filtering efficiencies reach- 
ing above 80% and 60% for Gr|Co and Gr|Ni interfaces, 
respectively. This work provides a unique analysis of the 
spin-polarized transport properties of a Gr|FM interface 
while considering the materials specific interactions in a 
non- equilibrium setting. 

The structural relaxations and band structure cal- 
culations were performed using DFT, the local den- 
sity approximation^ 1 for exchange-correlation poten- 
tials, the projector augmented wave method, and a plane 
wave basis with a cutoff energy of 400 eV as imple- 
mented in the Vienna ab-initio simulation package|12j. 
[jja) shows a diagram of the interface, where the left lead 




FIG. 2: Electronic structure of Co(lll)- and Ni(lll)- 
sandwiched graphene (left lead of system in[TJa)). The left 
and right panels show the majority spin and minority spin 
bands for Co and Ni respectively. The Fermi level is indi- 
cated by the horizontal dashed line. The black up-triangles 
(blue down-triangles) show the p z character of the top-site 
(hollow-site) C, while the red diamonds present the d z 2 char- 
acter of the Co or Ni atoms located directly above and below 
the top-site C. 



consists of FM-sandwiched graphene and the right lead 
is pure graphene. Transport occurs in the y-direction 
(as indicated by the red arrow) and periodic boundary 
conditions are assumed in the plane perpendicular to cur- 
rent, i.e., the system extends infinitely along the x- and 
z-directions. The most stable (minimal-energy) config- 
uration for graphene on a Co(lll) or Ni(lll) substrate 
was used [5] and corresponds to the top-site C located di- 
rectly above the FM atom and the hollow-site C sitting at 
the hollow site (see[ljb) and[ljc)). Seven FM layers are 
used to separate the graphene sheets in adjacent super- 
cells along the z-direction. The atomic structure of the 
graphene-FM contact was obtained by fixing the in-plane 
lattice constant to graphene's value of 2.46 A and relax- 
ing the atoms in the supercell until the net forces acting 
on the atoms were below 0.01 eV/A. The supercell box 
height, in the direction perpendicular to graphene, was 
varied after each relaxation in order to find the optimal 
height through total energy minimization. A fc-mesh of 
21x21x3 was adopted to sample the BZ for structural 
relaxations and total energy calculations. These param- 
eters provide the optimal graphene-FM distance do equal 
to 2.17 A and 2.13 A for Gr|Co(lll) and Gr|Ni(lll) re- 
spectively 

The spin-dependent band structure of the left lead is 
shown in [2j The green lines correspond to the bands of 
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the hybrid graphene-FM system. To locate the states 
originating from graphene, the weight of the C(p z ) or- 
bitals of all bands is superimposed (black up-triangle: 
top-site C, blue down-triangle: hollow-site C). There are 
three distinct features found from the figure: (i) At the 
K point, the graphene bands no longer show a linear dis- 
persion relation and now exhibit a band gap opening of rs 
0.4-0.5 eV, similar to what was shown for a single FM sur- 
face in contact with graphene [5] . The majority spin band 
gap resides at a lower energy than the minority spin band 
gap, and they do not overlap. This indicates that at spe- 
cific energies, one finds electrons of only a single spin type 
in the graphene. (ii) Superimposing the d z i character of 
the FM atoms located in the first layer above and below 
the graphene (plotted as red squares) shows that the FM 
interacts only via the top-site C. This can be seen from 
the fact that the red squares only overlap with the black 
up-triangles. Note that only the d z 2 orbitals showed any 
significant mixing with the graphene, in accordance with 
previous work [5]. (iii) It is also clear that the graphene 
conduction band (CB) and valence band (VB) result sep- 
arately from the top-site C and the hollow-site C, respec- 
tively. This is in contrast to pure graphene, where both 
C atoms in the primitive cell equally contribute to the 
CB and VB. The C-site dependence on the graphene CB 
and VB can be understood from the graphene-FM inter- 
action (mentioned in (ii)) which breaks the sub-lattice 
symmetry between the top-site C and the hollow-site C. 

Given the fully relaxed atomic structures and the well 
understood electronic states of the left lead, we em- 
ployed our state-of-the-art ab initio transport package, 
named MatDCal [T3], to compute the spin-dependent 
transport properties of the whole junction. MatDCal 
uses non-equilibrium Green's functions (NEGF) com- 
bined with DFT for open systems in a two-probe geome- 
try under finite bias, where the leads extend to ±oo. An 
optimized double-^ polarized atomic orbital set was built 
for each atomic species. The local density approximation 
[TT] and norm-conserving non-local pseudopotentials [H] 
were used and fc-point convergence tests, including the 
high symmetry points T, K and M, were performed for 
all calculations. 

The spin-dependent transmission coefficient of the mi- 
nority states (Tmin) in the 2D BZ (in the plane of k x 
and k z which are _L to the transport direction, i.e., y 
direction) for Gr(zigzag)|Ni is plotted in|3^a). The no- 
tation graphene(zigzag) or graphene(armchair) indicates 
the graphene is oriented such that transport occurs along 
the zigzag or armchair direction, respectively. Tmin 
shows a sharp peak at k x — 0, with near-zero values away 
from this point. To quantify this unusual dependence on 
k, we measure the width (Sk x ) of the Tmin peak. Sk x is 
defined as the half- max width of T, as depicted in[3ja). 
[3jb) presents the calculated spin-polarized 8k x as a func- 
tion of energy (E) for Gr(zigzag)|Ni, which clearly shows 
a nearly linear behavior. The spread in 8k x can be traced 
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FIG. 3: (a) Tmin versus kx at E = 0.1 eV for Gr (zigzag) |Ni. 

(b) Spin-polarized Sk x versus E for the same system as in 
(a). Sk x is defined as the width at half-max of T. The in- 
set illustrates the linear E-k dispersion bands of graphene. 

(c) Tmin versus k±_ at E = 0.1 eV for Gr(armchair)jCo. 

(d) Left panels: Top view in real space of Gr(zigzag)|FM 
and Gr (armchair) | FM. Right panels: Primitive BZ of pure 
graphene (zigzag) (top) and graphene(armchair) (bottom). 
The dashed arrows indicate the transport direction as well 
as delimit the zone edges along k x for our calculations, due 
to the use of a non-primitive supercell (the supercell width 
along x is shown in the left insets). The black curved arrows 
illustrate the folding of the K points. 



back to the conical states of pure graphene, which also 
exhibit a linear E-k relationship as shown in the inset of 
[3](t)) . For a Gr|FM interface, the pure graphene forming 
the right lead of our system only has electronic states 
at the K points (for E s» Ep). Hence, all the incom- 
ing carriers originating from the FM-covered graphene 
(i.e., the source or drain electrode) are blocked by the 
pure graphene except those with states at the K points. 
[3jc) presents Tmin versus for Gr(armchair)|Co. The 
armchair-oriented graphene also shows sharp peaks along 
k x , but are shifted in comparison to the zigzag-oriented 
graphene to k x = ±2/3. The positions of the peaks along 
k x can be explained, for both graphene orientations, by 
considering the effect of band folding on the graphene 
states. In[3]jd), we show the primitive BZ of graphene, 
where the dashed arrows indicate the transport direction 
in addition to delimiting the k x zone edges due to our 
non-primitive supercell. The K points located outside the 
k x zone edges will be folded inward, as illustrated with 
the black curved arrows. For graphene (zigzag), the con- 
ducting K points all appear at k x = leading to a single 
peak in Tmin , while the K points for graphene(armchair) 
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are positioned at k x = ±2/3, as seen in[3]ja) and^c). 

The fc-averaged spin-dependent transmission coef- 
ficients (T) versus energy E at equilibrium (zero 
bias) for Gr(zigzag)|Co(Ni) are presented in |4](a) 
and Qb), which were calculated using T„(E) = 
J BZ T a (E,k±) dk±, where Abz is the area of the 
BZ. For E smaller than s» Ep (Ep is set to zero), it 
is found that both Tmaj and Tmin vary roughly lin- 
early, similar to pure graphene. However, the T values 
for Gr|FM are roughly 50% smaller than those of pure 
graphene, a result of the band gap opening which reduces 
the band velocity. When considering E > Ep, one notice- 
able feature for both Gr|Co and Gr|Ni is the small Tmaj 
value between E « 0-0.4 eV. This energy range corre- 
sponds to the majority spin band gap of the graphene 
states (see [2]), thus explaining the decrease in Tmaj- A 
similar effect occurs for Gr|Ni but with Tmin (instead 
of Tmaj) between 0.43-0.80 eV, an energy range which 
overlaps with the minority spin band gap. |4|e) presents 
the spin transmission ratio (7) defined as T a /T s , where a 
is the spin component with the larger T value and a is the 
opposite spin. 7 is defined positive (negative) when a is 
the majority (minority) spin. Gr|Co and Gr|Ni both yield 
large dips near — 0.1 eV and 0.1 eV with 7 approaching 
18 and 25 respectively, each coinciding with the small- 
est T M aj value. Interestingly, for Gr|Ni, 7 shifts from 
minority spin dominated at 0.1 eV to majority spin dom- 
inated at 0.5 eV. This crossover behavior is attributed to 
the end of the majority spin band gap and the beginning 
of the minority spin band gap both located near 0.4 eV. 

The non-equilibrium calculations reveal the spin- 
polarized current (J (T )-voltage (V) characteristics of the 
Gr|Co and Gr|Ni interfaces, as plotted in [4^c) and (d). 
I (j is obtained from 

I ° = \\ T a {E)[f L {E,fi L )-f R (E,fi R )]dE, (1) 

where e is the electron charge, h is Planck's constant 
and f(E,/j,) is the Fermi-Dirac distribution. An applied 
bias V varies the left and right chemical potentials as 
Mi = Ef an d \ir = Ep + \e\V, where V = Vl — Vr. 
For Gr|Co, it was found that \Imin\ > I-Ta/ajI for all V 
values in the bias voltage window of interest. Whereas 
in the case of Gr|Ni, |/a//at| > \Imaj\ is only observed 
for positive V while the opposite result is found for neg- 
ative V (although less pronounced). To illustrate this, 
we plot the spin efficiency (17) for the Gr|FM interfaces 
in[4^f). 77 is calculated using \I a — I s \j\l a + I s \, where 
77 is defined positive (negative) when |/ma,/| > \Imin\ 
{\Imin\ > \Imaj\)- Gr|Co and Gr|Ni achieve maxi- 
mal spin efficiencies above 80% (at —0.2 V) and 60% 
(at 0.4 V) respectively, representing the percentage of net 
spin-polarized current. ?y — > for negative V in the case 
of Gr|Ni. This occurs because the integration window of 
E for calculating I a (ranging from (Ep + |e|V) — > Ep, 
valid for V < at zero temperature) is located below 
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FIG. 4: (a)-(b): Spin-polarized T versus E at V = for 
Gr|Co (a) and Gr|Ni (b) (majority spin: solid blue line, 
minority spin: dashed red line). The Fermi level is set to 
£7 = 0. (c)-(d): Spin-polarized / versus V for Gr|Co (c) and 
Gr|Ni (d) (majority spin: black circles, minority spin: green 
squares), (e) Transmission ratio 7 = T a /Ta (where a is the 
spin component with the largest T, and a is the opposite 
spin) versus E for Gr|Co (cyan full line) and Gr|Ni (purple 
dashed line). Note that T a fT„ is defined positive (negative) 
when T MAJ > T MIN (T MIN > Tmaj)- (f) Spin efficiency 
rj = \Ia — Ig\/\I a + Ig\ versus V for Gr|Co (cyan circles) 
and GrjNi (purple squares). Spin efficiency is defined positive 
(negative) when \Imaj\ > \Imin\ (\Imin > Imaj\)- 



both majority and minority spin band gaps (see [2]). This 
results in near-equal currents from both spin types. For 
Gr|Co, rj saturates slowly to zero with decreasing V, be- 
cause when the integration range extends to the bottom 
of the majority spin band gap, located near — 0.15eV, 
Tmaj remains small due to the vanishing density of 
states at the Dirac point in the pure graphene (which 
is pinned at the lower boundary of the integration win- 
dow). Hence, V must decrease beyond —0.15 V to obtain 
77 — > 0. The peak at V = results from Imaj and Imin 
both vanishing in pure graphene. The results shown in 
[4] have considered only graphene (zigzag). The transport 
properties of graphene(armchair) are found to be quali- 
tatively similar to graphene (zigzag). 

In summary, our non-equilibrium first principles trans- 
port calculations showed that Gr|Co(lll) and Gr|Ni(lll) 
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interfaces exhibit large spin injection values reaching 80% 
and 60% respectively. This effect originates from the 
graphene-FM hybridization which leads to an opening of 
the conical Dirac bands resulting in spin-dependent en- 
ergy gaps of the graphene states. Thus, in the ballistic 
regime, one can (theoretically) exploit the materials spe- 
cific bonding between graphene and the FM to achieve 
very efficient spin filtering. However, it remains to be 
shown whether these spin properties are robust in the 
inevitable presence of random interface disorder. 
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